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Abstract 

We compute the dynamics of excitation and two-body correlation for two-level "pseudoatoms" 
in a ID lattice. We adopt a simplified model where pair excitation within a finite range is perfectly 
blocked. Each superatom is initially in the ground state, and then subjected to an external driving 
laser with Rabi frequency satisfying a Poissonian distribution, mimicking the scenario as in Rydberg 
gases. We find that two-body quantum correlation drops very fast with the distance between 
pseudoatoms. However, the total correlation decays slowly even at large distance. Our results may 
be useful to the understanding of Rydberg gases in the strong blockade regime. 

PACS numbers: 03.67.Mn, 32.80.Rm 
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Recently there have been many experimental efforts in investigating Rydberg gases 



mm, 



spurred largely by fast developments in laser cooling and trapping. In such 
systems, the dipolar interaction between two nearby atoms will shift the pair excitation out 
of resonance with the driving laser. Local excitation is then greatly suppressed, showing 
the well-known phenomenon of dipole blockade. It is an important mechanism responsible 
for the explanation of many-body effects such as spectral line broadening and the 

expected sub-Poissonian atom counting statistics ^, Q, . It is also proposed to utilize dipole 



blockade to implement quantum gates useful in quantum computation and information 10|, 



HI]. To date, there have been many theoretical proposals and numerical simulations in 



exploring the excitation dynamics and coherence properties such as pair correlation and 
number correlation, which demonstrate the important role that dipole blockade played in 
these systems. 

Previous investigations of Rydberg gas have focused on the mean number of excitations 
and how the excitations are correlated in a gas. However, there hasn't been a study of 
genuine quantum correlation. This paper will discuss other useful information that can 
be extracted from the wave function, especially the two-body quantum correlation and 
total correlation (including classical correlation). As we will see later, two-body correlation 
is another important quantity for characterizing such systems. To investigate correlation 
properties, we are mainly interested in the following two questions. First, in such systems, 
how does correlation behave as the excitation begins to build up? Second, which correlation 
is more important? It has been shown that a mean field approach is not appropriate in 
dealing with correlation and a full quantum calculation is necessary 12|. However, the 
direct simulation of correlation is still absent due to numerical complexities in real Rydberg 
gases. 

We first consider the scenario as in current experiments where Rydberg atoms are confined 
in a volume V and driven by a laser. In the spirit of the dipole blockade model, we can divide 
V into equally sized regions V = IJ^ f a; with each region approximated as a sphere of blockade 
radius Rb. We follow the usual notation and call the ensemble of atoms in each region Vk a 



■'superatom" 



13|. For van der Walls interaction and a narrow band laser, is determined 



from Cq/R^ = Q, where Cq is the van der Waals coefficient in the atomic unit and Q is the 
Rabi frequency of the driving laser. Rydberg atoms are randomly distributed in V so that 
the number of atoms in each Vk is in general different. For example, in the homogeneous 
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case, the atom numbers in each superatom follow a Poissonian distribution with the same 
mean atom number. Such a picture is useful in explaining certain experimental results, e.g., 
the population dynamics under a driving laser. In this case, each superatom evolves with 
different Rabi frequencies. This leads to a fast and almost linear increase in the initial stage 
and saturation in the long time. However, such a model is not so interesting in showing 
the two-body correlation properties because correlation between different superatoms is 
apparently 0. Therefore, a proper modification to the superatom model is necessary in 
order to give nontrivial information on two-body correlation. 

In this paper, we discuss an elongated Rydberg gas whose transverse size is smaller than 
Rf,. This allows us to treat the Rydberg gas as a quasi ID system. The advantage is two 
fold. On the one hand, the edge effect in numerical simulation is smaller and the results 
converge relatively faster than those in higher dimensions with the same atom number. On 
the other hand, it permits an easier readout of the quantum state for ensemble of atoms 
in any region along the longitudinal direction. Similar to the 3D case, we divide such a 
quasi ID gas into a collection of superatoms aligned along the longitudinal direction. We 
investigate the two-body correlation by further dividing each superatom (vk) into smaller 
subregions labelled as Wk^a, i-e., Vk = [Ja=i'^k,a- The center of w^^a is denoted as r^a. The 
number of partitions in each Vk, labelled as N^, is assumed the same, Vfc. For the interest of 
this paper, we term each subregion Wk,a as a "pseudoatom" located at position a- in Fig. 
[H we give a schematic view of partitioning 3 superatoms into 9 pseudoatoms. The subscripts 
k and a can be further combined into a single one by relabelling them as in a spin chain. 
Interaction between two pseudoatoms on site j and k are assumed to be a van der Waals 
interaction with the distance \rj — r^l. This approximation neglects the position variance 
of excitation in each pseudoatom. By partitioning each superatom into more pseudoatoms, 
the error can be reduced. Different from the case of nearest neighboring superatoms (with a 
distance 2Rf,) where the interaction is usually ignored, we can see that nearest neighboring 
pseudoatoms (with a distance of a fraction of Rb) interact much more strongly so all of the 
pseudoatoms are correlated analogous to a spin chain. 

Due to the R~^ dependence of interaction, the transition from blockade to no blockade 
has a narrow width of R. This gives us the motivation to assume a perfect blockade with 
pseudoatoms within a certain distance and completely no blockade outside of it. From the 
many body point of view, the introduction of this approximation allows us to concentrate 
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FIG. 1: A schematic view of partitioning 3 superatoms into 9 pseudoatoms. Each Vk{k = 1,2,3) 
represents a superatom which is further divided into 3 pseudoatoms w^^aict = 1, 2, 3). 

on the much simpler short range interaction rather than the long range interaction. It not 
only lets us sec the effect due to pure blockade, but also greatly reduces the needed basis 
which makes the numerical simulation in a much larger system possible. The validity of 
this approximation will be discussed in the later part of this paper. Because of this type of 
partition for a superatom, the maximal number of pseudoatoms that can be blocked is just 
the maximal distance between two pseudoatoms j and k which have an interaction Vj^k much 
larger than the Rabi frequency Q, i.e., maxj^k{j — k\Vj,k/^ > v}- Besides the vagueness of 
this definition, in practice we choose rj = 6. A little algebra shows that the maximal number 
of pseudoatoms that can be blocked is given by maxfc=o,i,2...{^|^ < O.STN^}- The first few 
values are given here: for A^^ = 2, 3, 4, 5, 6, the maximal number of pseudoatoms that can 
be blocked is found to be 0, 1, 1, 1, 2, respectively. So the simplest nontrivial situation 
starts from Ny^ — 3, perfect blockade only between two nearest neighboring pseudoatoms. 
In this paper, we will limit our discussion to the situation where perfect blockade only exists 
between two nearest neighboring pseudoatoms. 

Because of our approximation, the Hamiltonian can be written as H — 
J2k=o JkP^^~^^x^P^^\ with the projector 

p(.) ^ -Q (1) 
q=k±\ 

where is the total number of pseudoatoms. a'^^'> are Pauli matrices for the fc-th pseu- 
doatom. — ^JXk/\ is the scaled Rabi frequency for the k-ih pseudoatom. We are 
assuming a resonant laser field with a constant intensity. is a random variable satisfying 
a Poissonian distribution with mean value A. In this case, A denotes the average number 
of Rydberg atoms in a pseudoatom. More exphcitly, A = N^/Nyj with Nj. the number of 
Rydberg atoms in a superatom which can be determined from the Rydberg gas density 
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n, the blockade radius Rb and the area in the transverse direction A through the relation 
Nr = 2nRhA. The underlying motivation for the assumed form of Jk is the number fluctu- 
ation and the collective excitation in each pseudoatom 6|. A given set of {Jq, Jn-i} is 
said to define a configuration. Our numerical results are obtained by averaging over many 
configurations. We note that the Hamiltonian formally describes a nonphysical multi-body 
interaction, which is of course an effective Hamiltonian constrained by the perfect block- 
ade requirement. All the pseudoatoms are initially in the ground state \g)^^ , and then 
subjected to a resonant external laser field. The subsequent dynamics constrained by the 
perfect blockade are what we are after in this paper. 

For nearest neighboring blockade, the maximally possible excitation for all pseudoatoms 
is found to be where [■] denotes the integer part. As a consequence, the number of 

restricted basis, Nb, is significantly reduced from that of the full basis set. To compute the 
wave function, we expand it in the restricted basis labelled as |/ip), i.e., 

^W = 5Zcp(t)|^p), (2) 
p=i 

the fraction of excitation is found to be 

^ Nb N-l 

Pe.it) = ^5^|cp(t)|'(/xp| J2 Ie)fc(e||/ip). (3) 

p=l k=0 

In Fig. [21 we show the numerical results of excitation fraction as function of time for 
pseudoatom number = 16. The results are obtained by using wrap boundary condition 
and averaging over 1000 different configurations of {Jk}- We can see that the overall trends 
of the curves for different A are similar. In all cases, Pg^ first increases almost linearly 
and overshoots to a maximal value, then oscillates and saturates to about 26%, close to 
the expected value (25%). It comes from the fact that for two pseudoatoms with perfect 
blockade, the wave function should behave like cos{Jt)\gg) + sm{Jt){\eg) + \ge))/\/2, so the 
total time-averaged excitation for one pair of the pseudoatoms is 50%, i.e., 25% for each 
pseudoatom. It's also obvious that the larger A, the less fluctuation in Jk, so that we can 
observe more oscillations. For A = oo, when there is no fluctuation at all, the oscillation 
period can be estimated from analyzing the spectrum of the Hamiltonian. Labelling the 
eigenvalues and the corresponding eigenstates as Ea and \Ea) (a = 1,2,..., Nb) respectively, 
we find that there exists a nearly periodic structure in the plot of Ki'it = 0)\Ea)\'^ versus 
Ea- In this case, Pex does not saturate to a fixed value for A^ = 16 pseudoatoms. 
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FIG. 2: Excitation fraction P^x as function of time for different A. The inset shows enlarged 
version of initial dynamics. 

The excitation fraction of a physical Rydberg gas can be quahtatively obtained by mean 
field approach j^, which is essentially a reduction to single-atom picture. However, the 
two- atom correlation, where the off-diagonal terms of the density matrix are important, 
cannot be simply accounted for from the mean field approach. In the following, we will 
focus on correlation properties of reduced two-pseudoatom subsystem from a full quantum 
mechanical calculation. 

The first quantity we will look into is the pair correlation. Pee, at saturation time, which 
is defined as the ratio of the probability of both pseudoatoms excited divided by the square 
of probability of single-pseudoatom excitation 12]. In Fig. [3l w e show Pee as function of 
distance for different A. Compared with the results in Ref. 12] (see Fig. 2) which uses 
a continuous model, our results captures the main physics in a much simpler way. We 
can see that for both the discrete and continuous model, P^e exhibits similar behaviors: it 
is negligible within a distance and there is a sharp increase after that; at large distance, it 
saturates to 1, corresponding to no correlation. For all three A in our calculation, the next to 
nearest neighbor always has the largest value. This can be interpreted as follows. Assuming 
the index of the first pseudoatom to be 0, the separation d is simply the label of the other 
pseudoatom. Pee for (i = 1 is because pseudoatom 1 cannot be excited once pseudoatom 
is excited. For d = 2, P^e is well above 1. This is because if pseudoatom is in the excited 
state, then pseudoatom 1 must be in the ground state. Thus pseudoatom 2 is more likely 
to be in the excited state. Similarly, pseudoatom 3 will be less likely in the excited state. 
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resulting in a value lower than 1. The fast approach to 1 at large distance means that such 
pair correlation is only short ranged. We do not present the result for A = oo due to the 
non- negligible oscillation even at long enough time. 
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FIG. 3: Pair correlation Pge as function of distance for different A. Solid lines are guides to the 
eye. 

While the above pair correlation only relies on the diagonal elements of the reduced 
two-body density matrix, the genuine quantum correlation must be computed from the full 
reduced two-body density matrix. We use entanglement of formation (EOF) as our measure 
for quantum correlation [l^. It is related to the so-called concurrence, C = max(0, a/Ai — 
— "s/As — y/Xi) , where Aj are the eigenvalues of a Hermitian matrix, R = ^{py ® 
ay)p*{ay ® (yy)\fp-, in descending order. EOF [e) is then given by e = /i((l + vT^^C^)/2), 
where h{x) = — xlog2(a;) — (1 — x)log2(l — x). e gives for separable states and 1 for 
maximally entangled states. 

In Fig. m we show EOF for pseudoatom pairs (0, 1) for different A. We find that 
shows similar pattern as excitation fraction. During the same time scale, they approach a 
maximal value and then decay slowly. This reflects two simple facts that entanglement is 
originated from excitation and it is superposed coherently in the initial stage. The behaviors 
of entanglement for other pairs are shown in Fig. O We can see they again behave in a 
similar way, although the first peaks are offset in time. 

To investigate how entanglement varies with the distance between pseudoatoms, we use 
the first peak of EOF for different pairs to characterize its dependence on the distance. Our 
results are shown in Fig. [61 We note that EOF drops almost exponentially as the distance 
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increases for all A. Thus we emphasize an important point here: only nearest neighboring 
pseudoatoms have relatively large entanglement. Entanglement of other pseudoatom pairs 
are negligible during the full time scale we investigate. This is consistent with our intuitive 
understanding that there is no long-range order in such systems. 
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FIG. 4: EOF e as function of time for different A. The inset shows enlarged version of initial 
dynamics. 
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FIG. 5: e as function of time for different pair with A = 3. 



We can see that quantum correlation is significant only for nearest neighboring pseu- 
doatoms. However, this is not the case for the total correlation (including classical correla- 
tion). To compute the total correlation, we adopt the measure as suggested by D. L. Zhou 
and L. You [15|]. For density matrix p(i2) of pseudoatoms 1 and 2, it is given by 
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M, 



-Tr|p(i2) -P(i) ® P(2)|, 



(4) 
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FIG. 6: The first peak of EOF as function of distance for different pair. Solid lines are guides to 
the eye. 



where p(^k) is the reduced density matrix of the k-th pseudoatom. Different from Ref. 15| . 
here the prefactor 2/3 is chosen so that this measure gives 1 for a maximally entangled state. 

In Fig. [71 we show the total correlation as function of time for different A. We again 
see that they show similar pattern as EOF. To investigate the distance dependence of total 
correlation, we show the first peak value of total correlation in Fig. [HI We can see that 
total correlation does not drop like a single exponential function with distance as that in 
EOF. Rather interestingly, they even do not decrease monotonically as increasing distance. 
A simulation with = 20 atoms shows similar behaviors with distance. We find that 
they give almost identical results for several initial distances. However, close to the middle 
point of the lattice (distance N/2), those values of = 20 have correspondingly lower 
values than A^ = 16. Therefore, such non-monotonicity will become negligible for A^ oo, 
i.e., a finite size effect. We conclude that although quantum correlation is not important 
in most cases, the total correlation cannot be simply ignored, i.e., there is non-negligible 
classical correlation in such systems. It is necessary to include not just nearest neighboring 
pseudoatoms when considering correlation-related properties. 

The extension of our results to higher dimensions seems difficult due to the computational 
limitation. For 2D configurations and nearest neighboring blockade, we can only simulate 
up to 5 X 5 square lattice. Without presenting more figures, we briefly discuss our results. 
Our major results of ID still hold in 2D configurations. We find that EOF shows similar 
patterns but drops even faster as distance increases. This is possibly due to the fact that in 
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FIG. 7: Total correlation of pair (0, 1) as function of time for different A. The inset shows enlarged 
version of initial dynamics. 
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FIG. 8: The first peak value of total correlation as function of distance. Marked curves are for 
= 16 with corresponding A shown in the legend. Unmarked curves are for = 20 with the 
lower (dashed) and upper (dotted) curve corresponding to A = 3 and A = 48, respectively. 

2D, as more neighboring pseudoatoms are involved, the fluctuation is more intensive than 
that in ID. In addition, the corresponding pair entanglement is found to be smaller than 
that in ID. Thus we conjecture that the pair entanglement is even smaller in a 3D cubic 
lattice. 

Before concluding, we want to justify the validity of our simplified Hamiltonian. We will 
show that the fast decay of entanglement on distance is not an artifact of the assumed short 
range interaction. To confirm this point, we carry out simulations using the long range 
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interaction, 



D 



1 + ai^-) 1 + a?) 



b - k\ 



(5) 



fc i<fc 

where Jk takes the same form as in the Hamiltonian H . The dimensionless parameter D 
quantifies the interaction strength between two nearest neighboring pseudoatoms. Because 
we always use the distance between two nearest neighboring pseudoatoms as the length 
scale, D also depends on A^^, i.e., the decrease in the distance is equivalent to the increase 
in the interaction strength. We carry out similar calculations by using the full basis of H' . 
Selected results for A^^ = 5 with nearest neighboring blockade are shown in Figs. [9] and [TOl 
for Smax and Mc^^ax as function of distance, respectively. We find that the main features 
by using H' are qualitatively reproduced by the simplified Hamiltonian H , which leads us 
to conclude that the behavior of entanglement on distance is not an artifact of short range 
interaction in the simplified Hamiltonian H. And for the total correlation, the overall shape 
of the curves remains the same. More importantly, the total correlation does not drop as 
fast as entanglement with distance, which provides another evidence that our simplified 
model captures the main physics. We emphasize that this calculation not only justifies our 
simplified model, but also means that those effects in Figs. El [HI E], and [lO] are mostly due 
to blockade, which shows that the long range interaction is actually not so important to 
correlation properties. 
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FIG. 9: £max function of distance for Ny^ — 5 with nearest neighboring blockade. 

In conclusion, we have performed full quantum calculations on correlation properties for 
pseudoatoms in a ID lattice structure with perfect blockade. By comparing the results from a 
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FIG. 10: Mc max 3-8 function of distance for N^^ — 5 with nearest neighboring blockade. 

reduced basis and a full basis calculation, we justify the validity of the simplified Hamiltonian 
for perfect blockade among pseudoatoms. This agreement also means that single- and two- 
particle correlation properties (e.g., average excitation, correlation) are determined mostly 
by the pure blockade effect. From numerical simulation, we find that there are both quantum 
and classical correlations accompanying the building up of atomic excitation. Our results 
show that two-body entanglement is only important for nearest neighboring pseudoatoms 
and it drops exponentially fast with the distance between them even when there is no 
fluctuation in the system. However, the total correlation decays much more slowly with 
distance, showing the system in this paper is mostly classically correlated. As a simple 
extension to higher dimension, we compute our model system in a two dimensional 5x5 
square lattice. We find 2D results agree qualitatively with those of ID. From the theoretical 
point of view, our findings imply that a better description of Rydberg gas beyond mean 
field or superatom picture should also at least take classical correlation into consideration. 
We hope that our study can be helpful to the understanding of Rydberg gases in the strong 
blockade regime. 

We would like to thank J. V. Hernandez for some helpful discussions. This work is 
supported by the NSF under grant no. 0355039. 
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